Automatic Development of Global Phase Diagrams for Binary Systems in Pressure- Temperature Space

In presenting this thesis in partial fulfillment of the requirements for a Postgraduate degree from the University of Saskatchewan, I agree that the Libraries of this University may make it freely available for inspection. I further agree that permission for copying of this thesis in any manner, in whole or in part, for scholarly purposes may be granted by the professor who supervised my thesis work or, in his absence, by the Head of the Department or the Dean of the College in which my thesis work was done. It is understood that any copying, publication, or use of this thesis or parts thereof for financial gain shall not be allowed without my written permission. It is also understood that due recognition shall be given to me and to the University of Saskatchewan in any scholarly use which may be made of any material in my thesis. ABSTRACT Global phase diagrams of binary systems in pressure-temperature (PT) space are very useful. In this project the techniques to automatically develop global phase diagrams in PT space were created. The codes to compute different components of a global phase diagram in PT space were developed. These codes were then successfully incorporated into a single functional program. To generate the binary PT phase diagram, the overall composition was varied from pure component 2, the least volatile component (LVC) to pure component 1, the most volatile component (MVC). The step size for changing mole fraction was varied in the calculation of different parts of a global phase diagram. When the points near the joining points between different parts were computed, the step size was set to a rather small value. The step size was then increased to twice of the last value for each subsequent point computed. When the MVC mole fraction was approaching one, the step size was set to a small value to obtain enough points needed to minimize the chances of missing important phenomena. The techniques to set initial guesses for evaluation of different components of a global phase diagram were discussed. The code performance, including the number of iterations for different convergence criteria and the sensitivity of the algorithm were presented. Using the code developed, phase diagrams of type I, type II, type III and type V were generated using representative binary systems from the petroleum processing field. iii The boundary states between different types of phase behaviour were also …